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Abstract 

We review past work using broad emission lines as virial estimators of 
black hole masses in quasars. Basically one requires estimates of the emit- 
ting region radius and virial velocity dispersion to obtain black hole masses. 
The three major ways to estimate the broad-line emitting region (BLR) ra- 
dius involve: (1) direct reverberation mapping, (2) derivation of BLR radius 
for larger samples using the radius-luminosity correlation derived from rever- 
beration measures, and (3) estimates of BLR radius using the definition of 
the ionization parameter solved for BLR radius (photoionization method). 
At low redshift (z £ 0.7) FWHM H/3 serves as the most widely used estima- 
tor of virial velocity dispersion. FWHM H/3 can provide estimates for tens 
of thousands of quasars out to z m 3.8 (IR spectroscopy beyond z « 1). A 
new photoionization method also shows promise for providing many reason- 
able estimates of BLR radius via high S/N IR spectroscopy of the UV region 
1300 - 2000 A. FWHM MgllA2800 can serve as a surrogate for FWHM H/3 
in the range 0.4 ^ z £ 6.5 while CivA1549 is affected by broadening due to 
non-virial motions and best avoided (i.e. there is no clear conversion factor 
between FWHM H/3 and FWHM CivA1549). Most quasars yield mass es- 
timates in the range 7 £ logM B H ~ 9.7. There is no strong evidence for 
values above 10.0 and there may be evidence for a turnover in the maximum 
black hole mass near z ~ 5. 
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1. Estimation of Black Hole Masses in Quasars: an Introduction 



The existence of supermassive compact objects in galactic nuclei was 
pr oposed to solve the energy problem raised by the discovery of radio galax- 



ies 



Hoyle and Fowled (119631 ) . They estimated that the mass of such an ob 



ject had to be of the order of 10 8 solar masses. iHoyle and Fowlerl (119631 ) 
showed that nuclear reactions were insufficient to provide the energy and that 
the energy source had to be g ravitational collapse (see a lso ISalpete r. 1964; 



Zel'Dovich and Novikovl . Il965l ) . The d iscovery of quasars ( iGreenstein and Schmidt 



19641 ) and of rapid optical variability ( ISmith and Hoffeitl . Il965l ) exacerbated 



the problem because, at their redshift distances, many quasars apparently 
involve sources emitting 10 3 times the luminosity of an L* galaxy within 
a volume much less than a parsec of diameter. Estimation of the masses 
for these black holes is therefore at the center of quasar astrophysics and 
cosmology. 

The black hole mass of quasars is a fundamental parameter that relates 
to the evolutionary stage of quasars and of the accretion processes occurring 
within them. An estimate of black hole mass (Mbh) allows one to assess 
the role of gravitational forces in the dynamics of the region surrounding the 
black hole. The power output of quasars is directly proportional to M B h- 
There is much debate over how the evolution of quasar energetic output 
might affect the development and structure of the host galaxy, as well as 
larger scale structure formation a nd, at very high redshifts, the re-ionization 



of the Universe (e.g., iFanl . |2010| ). 



Early estimates of Mbh were based on the apparent similarity of quasar 
spectra and their line widths. The size of the emitting region can be written 
as: 



^ ( _3£line_\ 3 /js 

where L\\ ne is the line (typically H/3) luminosity, /f is the filling factor 
and eii ne is the line emissivity. This expression was used to derive Mbh by 
considering the nucleus as a bound system and line broad ening due to gas 
cloud motions (IDibail . 119771 Il984t IWandel and Yahill . I1985I ). It is interesting 
to point out that an analogous argument (coupled with the relatively short 
timescale of observed continuum variations) was used immediately after the 
discovery of quasars to suggest a small emitting region size and large mass. 
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The similarity of quasar spectra implies a roughly constant i onization pa- 



rame t er U or a constant product U times the electron density n e ([Baldwin and Netzerl . 



19781 ; iDavidson and Netzerl . Il979l ) , where 



U = Ju ° h \ . (2) 
A-Kn e cr 2 

Here L u is the specific luminosity per unit frequency, h is the Planck 
constant, z/ the Rydberg frequency, c the speed of light, and r can be inter- 
preted as the distance between the central source of ionizing radiation and 
the line emitting region. So basically tblr is oc L 1 / 2 using rather order-of- 
magnitude considerations and the assumption of constant U pointing towards 
large masses increasing with source luminosity. 

2. The Virial Assumption 

The virial mass is defined as: 

M BH = f-^, (3) 

where / is a factor dependent on geometry of the emitting region, r 
is the distance of line emitting gas from the central black hole, Sv T is the 
line broadening due to virial motions, and G is the gravitational constant. 
All methods based on broad optical and UV lines (e.g. photoionization, 
reverberation) come down to estimating the broad line region (BLR) distance 
from the central continuum source (hereafter indicated as the BLR radius 
?"blr,) and combining it with a broad line width to derive Mbh- 

The term "Keplerian" is used to indicate the particular case of a gravita- 
tionally bound system where the virial theorem applies. Kepler's third law is 
appropriate when mass is mainly concentrated in a single body and the self- 
interaction gravitational potential of any remaining mass can be neglected; 
in the case of quasars, Mbh^> mass of the line emitting gas (i.e., not unlike 
the solar system if solar and planetary masses are compared). Therefore, the 
terms virial and Keplerian will be used synonymously in the quasar context. 

The fundamental question is whether the virial assumption is correct for 
the line being used. Evidence in support of vi rialization come from early ve- 



locity resolv ed reverberation mapping studies (lGaskelll . ll988uKoratkar and Gaskelll . 



19891 . Il99ll ) that excluded outflow as the broadening source in at least low- 



ionization lines. Circumstantial evidence for virial motion involves emission 
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Figure 1 : Test of virial broa dening on several emission lines of NGC 5548. Adapted from 
Pete rson and Wandell (Il999h . The line shows an unweighted least square fit. 
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l ine profiles of NLSyl-li ke sources which are relatively symmetric and smooth 
flMarziani et all l2003ah . 

If the virial assumption is valid then the velocity field is Keplerian and 
l ine broadening should anti - correlate with the time lag of dif ferent lines 
flPeterson and Wandell . Il999l : IPeterson et all liooi iKrolikl . bOQlh . Sufficient 
data to test the width-distance relation are available for NGC 5548. Fig- 
ure [T] shows a trend (r oc 8v^) consistent with dependence of time de- 
lay on line width. The delays of 6 emission lines (including the blend 
Olv]A1402+SilvA1397) are plotted versus FWHM measures. Data come 
from the International AGN watch carried out in 1989. An unweighted least 
squares fit yields a slope of a = 2.01 ±0.60 which is consistent with the value 
expected for a Keplerian velocity field. A weighted least squares fit that 
takes into account errors on both axes yields a slope 2.60 ± 0.60 indicating a 
strong increase in FWHM with decreasing distance. 

Equation [3] is deceptively simple. Every term is problematic: 

1. The geometric factor / allows us to convert an observed FWHM value 
into a virial velocity. / is poorly known because BLR geometry and 
effects of line-of-sight orientation are poorly known. 

2. The distance r is highly uncertain because the BLR is not spatially 
resolved in any source. 

3. The appropriateness of 5v (virial broadening) rests on the assumption 
that lines are Doppler broadened reflec ting motions of the emitting gas, 
and not by scattering (as considered by Mathis 1970l ; Kallman and Krolik 
19861 ; iGaskell and Goosmannl 120081 ) . To estimate 5v one must deter- 



mine which line or line component (if any) is broadened by the motions 
of virialized gas. Line profile shifts relative to the quasar rest frame 
warn against indiscriminate use of line profile widths to estimate virial 
broadening. 



3. The factor / 

Our uncertainty about the geometric parameter likely represents a road- 
block to estimation of black hole masses more accurate than a factor of 3 - 
5. If FWHM of a broad line is chosen as an estimator of Sv T and the velocity 
distribution is isotropic then the square of the velocity module is v 2 = 3v^ 
implying a geometric factor / w a/3/2. In Eq. [3] our lack of knowledge 
about the geometry and dynamics of the BLR is concentrated in /. In other 
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words, an understanding of the dynamics and geometry of the BLR is di- 
rectly related to the determination of M B h- Considerable recent discussion 
and modeling have been directed tow ards estimating f and how it might de 



pend on changing physical conditions ( ICollin et al.l . l2006t lOnken et all 12004 



Netzer and Marzianil . 120101 ) . Radiation press ure forces and orien tation effects 
make it possible that it is significantly (x3 in lCollin et al.ll2006l ) and system- 
atically different in Pop. A (FWHM H/3 < 4000 km s" 1 ) and Pop. B sources 
(FWHM H/3 ^ 4000 km s _1 ). Pop. A and B show striking di fferences in the 



H/3 a nd CivA1549 profiles and in several other properties (jSulentic et al. 
20071 ). Although their interpretation in terms of BLR structure is not yet 
clear, further discussion in £0 will enlighten the importance of separating 
quasars into these two populations. 

3.1. Dynamical effects 

Given the strong radiation field to which gas is exposed in quasars it 
seems unlikely that gravity is the only important force. The relative balance 
of gravitational and radiation forces in a source of fixed Mbh points toward 
a role for Eddington ratio. The ratio between radiative and gravitational 
acceleration in an optically thick, Compton-thin, medium can be written as 



^rad 



7 o Lho1 AT- 1 

' - Z T JV c,23 
Lit: • 



(4) 



N, 



c,23 



grav -"-^Edd 

where L bol and L Edd are the bolometric and Eddington luminosities 
is the hydrogen column density in units of 10 23 cm" 2 . If r a , > 1 
five a cceleration dominates (cf. Netzer and Marziani . 2010l ; Marziani et al. 
2010h . 



Very high column density clouds will be unaffected by radiative forces 
while lower column densities will be affected even at moderate Eddington 
ratios. A more refined dynamical treatment should consider the effects of 
radiative forces on a system of clouds moving in pressure balance with an 
external medium. 

This approach indicates that the effect of radiation pressure significantly 
alters the equilibrium of a system of clouds: the emissivity weighted radius 
will increase as orbits become more elliptical and hence the cloud will spend 
a longer time at larger distances from the central continuum source. The 
value of / will therefore depend on Eddington ratio. At the same time the 
line will become narrower as radiation pressure increases. Fig. [2] shows 
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two H/3 profiles computed for confined clouds moving in pressure equilib- 
rium with an external med ium assuming a spherical distribution of orbits 
(INetzer and Marzianil . I201CH ) . Eddington ratio that has been set to L/L-Edd= 
0.05 (at this low L/L^dd the effect of radiative forces is small) and to L/L~Edd= 
0.5 (appropriate for NLSyl-like Pop. A sources). The effect on FWHM and 
on the virial product tblrFWHM 2 is rather modest, and the derived / in- 
creases from ^0.76 to w 1.08 along with L/L Edd . How will these results affect 
single-epoch Mbh estimates in a large sample? The dependence of tblr on 
Eddington ratio will be a source of scatter in the tblr-L correlation. If this 
is taken into account we find an effect of less than a factor of 2 for Mbh 
estimates over 4de x in luminosity Netzer and Marziani . 2010l ). 

Previous work ( Marconi et al. , 2008 . 20091 ) accounted for radiation pres- 
sure effects by adding an Mbh term proportional to the ratio L/N c (where 
N c was supposed to have a log- normal distribution with average logiV c = 23). 
In this case corrected Mbh can be a factor of ev en £ 10 higher than values 
obtained using the virial formula (jGaskelll . Il996l ). The implied masses of the 
most luminous quasars that would have been ^ 10 10 M Q , p ossibly r eintro - 
ducing what we have called the "maximum mass problem" (iNetzerl . 120031 ): 
the most massive black holes at any redshift cannot exceed by a large factor 
the maximum mass of the (mostly quiescent) black holes at low z, aroun d 
5 • 10 9 M (the mass of the black hole in Messier 87; iMacchetto et al][l997j ). 
In ot her words, black hole s with masses above 5 • 10 9 M should be ext r emely 
rare (ISulentic et al.l 120061 and references therein) . The iMarconi et al.l (120091 ) 
results are probably incorrect because they neglect the likelihood that the 
column density of bound clouds is dependent on Eddington ratio and hence 
on luminosity for a given mass. 



3.2. Estimating f empirically: normalization 

The zone of gravitational influence for the black hole has been resolved 
in several near by galaxies. The most reliable methods involve observation 



19951 : iGreene et all l2010bl : iKuo et all 1201 if ) 



of n 2 masers (IMiyoshi et al. 
or, at a s econd stance, space-based long-slit sp e ctra of optical emission 



lines (e.g.. iFerrarese et al.l . Il996l ; IMacchetto et al.l . Il997l ; iBarth et all 12001 



Capetti et al.l . 120051 ). Both methods spatially resolve the velocity field and 



yield a radial velocity "cusp" indicative of distances where the velocity field is 
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Figure 2: Effects of radiation pressure on l ine profiles computed f or a s pherical distribution 
of clouds in pressure equilibrium following iNetzer and Marzianil (|2010l ). Model parameters 
are kept the same save for the Eddington ratio that has been increased by a factor 10, 
from L/LEdd= 0.05 to L/L^dd= 0.5 . The dot-dashed lines mark the rest- frame origin 
and the half maximum intensity level. 



governed by gravity of the black hole. Mass values obtained from gas dynam- 
ics are considered less reliable since the line emitting gas might be more af - 
fected by non-circular motions (e.g. , iBarth et al.l . l200ll ; IShapiro et al.l . 120061 ) . 
Black hole masses computed from resolved velocity fields c orrelate with 



stellar velocity dispers i on me asures in the bulges of host galaxies ( iFerrarese and Merritt 



2000l ; iGebhardt et al.l . 120001 ) . The Mbh - correlation obtained for veloc- 
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ity resolved galaxies has been used to estimate / which can then be used 
for computation of quasar BH masses obtained from reverberation mapping 
derivations of tblr- One basically overlaps the Mbh - ®± relations for non- 
active galaxies and AGN where a* has been estimated, usually from the IR 



Woo et al. 



2010 



Graham et al. 



Calcium tr i p let ( Onke n et all 12004 ICollin et all 120061 : iGiiltekin et all |2009; 



2011). 



The derived f a ~ 5.5 of lOnken et all ( 120041 where j a refers to the / 
value when the velocity dispersion is used as virial broadening estimator) 
was reputed a valid estimate but a recent study now suggests f a ~ 2.9 
( jGraham et alll201ll ). halving the Mbh values. Considering that the Mbh _ 
a* relation for nearby galaxies shows significant scatter, results are subject 
to substantial uncertainty. In addition, there are two fundamental caveats: 

1. Mbh - (J* c orrela tion of nearby galaxies could be a selection effect 



( IKormendyl . Il993l ): only the most massive black holes can be resolved 
for a given bulge mass. The relation between M B h and a* should be 
taken with special care in the lower Mbh range. This issue is farther 
discussed in §|9j 

2. It is unlikely that a single / value is valid. / appears to d epend on the 
width of the H/3 line or, almost equivalently, on its shape. ICollin et al. 
(120061 ) separ ated the reverb e ration sample in two populations: the Pop. 
A and B of ISulentic et all (j2000al ) and Pop. 1 and 2 according to a 
profile shape criterion. This difference reflects different geometries that 
should also lead to di fferent respo n ses to continuum changes. This most 
interesting result of ICollin et all (120061 ) should be taken into account 
in any eventual attempt to derive /. 



3.3. Orientation Effects 



Orbits of BLR clouds are unlikely to be oriented randomly resulting 
in a strong depen dence of M B h on orientation. The / value derived by 



Collin et all (120061 ) indicates that low-ionization line (LIL) emission in Pop. 



A sources may arise in a flattened configuration. Other lines of evidence sug- 
gesting a flattened gas distribu tion have been reported since the 1970s (they 
have been recently reviewed by iGaskelll l2009bl ) . Sources viewed at high incli- 
nation will then be measured with too small FWHM leading to systematic 
underestimation of Mbh- 
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Figure 3: Effects of orientation on the low-z sample of Marziani et al. ( 2003bl ). Blue 
outliers (large blue spots) are indicated along with other sources in our sample (filled 
circles) in the plane L/L^d vs. A/bh- The largest circle indicates the position of PKS 
0736+01, a radio- loud blue outlier. Arrows indicate displacements of blue outliers if we 
apply an orientation correction to their masses. 



Measures of both radio-quiet and radio-loud quasars likely suffer from 
inclination effects with the latter offering a rather straightforward way of 
estimating its amplitude. One can compare average FWHM H/3 values for 
core-dominated and lobe-dominated sources. If the radio jet is oriented per- 
pendicular to a flattened distribution of line emitting clouds then core- and 
lobe-dominated represent sources where the disk is viewed near and far from 
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face-on r espectively. Lobe-dominated sources are found to be broader by a 



factor 2 flMilev and Milled . Il979l : ISulentic et all . 120031 : IZamfir et all . |2008|). 



The assumption of an isotropic velocity component combined with a rota- 
tional one roughly envelops the dist ribution of radio-loud sou rces in the plane 
R vs FWHM H/3 parameter plane ( jWills and Brownd . Il986l ). In a small sub- 
sample of radio-loud sources with detected superluminal motion the apparent 
speed can be used to retrieve the angle 9 between the line-of-sight and the 
relativistically moving plasma, once an estimate of the Lorentz factor 7 is ob- 
tained in the fr amework of the synchrotron self-Compton model of X-ray and 
radio emission (jRokaki et al.l . l2003t ISulentic et al.l . 120031 ). A diagram of 9 vs. 
FWHM confirms a factor ^ 2 broadening for the sources observed at larger 
inclination. An effect of similar amplitude is expected for radio-quiet sources, 
although attempts have yet to yield a con vincing; orientation indicator (see 
however the recent work of iBorosonl 120111 ) . Considering that 9 may be in 
the range few degrees ^ 9 £ 45°, and distributed randomly, we can have an 
average < 9 >m 30°; this means that, ignoring orientation effects leads to an 
underestimate of the mass for almost all sources, up to a factor ~ 10 if the 
source is observed almost pole-on and if the BLR is strongly flattened. The 
effect is less (a factor ps 5) in the case of Pop. B sources where a vertical 
component of motion (p ossibly due to turbulence) is more strongly affect- 
ing; the line width (e.g., iGaskeliliooibh . Even if pole-on sources are rarest 
in a randomly-oriented sample, errors of this amplitude may dramatically 
increa se the scatter in large sam ples and influence inferences on Eddington 
ratio (I Jarvis and McLurd . 120061 ) . Fig. [3] shows the rare population of "blue 
ou tliers" i.e., Po p. A sources showing a large [Qui] AA4959,5007 blueshift 



(Za manov et al 



20021 ) and believed to be oriented almost pole-on (IBoroson . 



201ll ). If broad H/3 is emitted in a highly flattened configuration then they 
will appear systematically undermassive and hence super-Eddington. An 
orientation correction of 0.4 (in logarithm) to estimated masses of the blue 
outliers moves them in the sub-Eddington regime of LjL^A vs Mbh diagram. 



4. BLR radius from Reverberation Mapping 



Gaskelll . Il988h 



The last 25 years (the first recorded attempt is due to 
saw the rise of reverberation mapping based estimates of Mrh arid many re- 
views can be found on this technique (jPetersonl . Il993l ; iPeterson and Hornd . 
20061 ; iPetersonl . 120081 ) . The cross-correlation function bet ween continuum and 
emission-line light curves yields a time lag measure t (IGaskell and Sparkd . 
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Figure 4: Correlation between t ime delay in units o f day and luminosity at 5100 A in units 
of erg s" 1 . Original data from lBentz et al. ( 2009t h separate determinations for the same 
object are shown as independent points. The filled (blue) circles refer to NGC 5548; the 
dashed line is an unweighted least square fit. The filled line is a least-square fit weighted 
over the uncertainty in logt. 



1986T ). The time lag due to the light trave l time across the broad line emit- 
ting region yields an estimate of tblr= ct (ILyutyi and Cherepashchukl . Il972 



12 



Cherepashchuk and Lyutyil . Il973l ) . One can then derive the black hole mass 
using FWHM H/3 B c under the assumption of virialized motions, using Eq. El 
Here we will focus on the accuracy of such Mrh estimations. The measure of 
rm/R is based upon several assumptions flGaskell and Sparkd . Il986l ; iKrolikl . 



20011 ; iMarziani et all 120061 ). 



2. 



The continuum emitting region is much smaller than the line emitting 
region. This assumption might not be satisfied since the width of the 
continuum autocorrelation function is usually not much less than the 
measured delay. If the continuum is emitted by the accretion disk, the 
region emitt ing the optica l continuum might reach the inner BLR (see 
Table 1 of Gaskell . 20081 ) . However, the autocorrelation function of 
the ionizing continuum (usually not observed) might well be narrower 
than the one considered in most reverberation studies (see e.g.. lGaskelll . 



20081 . for a review). For H/3, the continuum might be monitored with 
broad band photometry, or one may measure the specific continuum 
flux at a suitable wavelength; 

The observed continuum and the Hi ionizing continuum are directly 
related. The latter assumption appears to be valid, since the monochro- 
matic luminosity at selected UV wavelengths has been shown to strongly 
correlat e with the BLR radius by the international AGN Watch c am- 
paig ns (Clavel et al. 1991 ; Peterson et al. 1991 ; Edelson et al. 19961; see 



also iMcLure and Jarvisl 120021 ; IVestergaardl l2002l and iGaskej \20m for 
an analysis of NGC 7469, the only object for which a significant time 
delay between the optical and UV continuum was found). 
The light travel time across the BLR is shorter than the dynamical 
response time t dyn (so that BLR structure does not change over the 
light travel time). This has been verified by the monitoring campaigns. 
No measurable dynamical effects due to radiation pressure are present. 
This is probably the case for relatively short time lapses (< 1 yr in 
nearby Seyfert 1 galaxies) i.e., a timescale shorter than tdyn ~ ^blr/ 
Av r ~ 4 yr in the case of NGC 5548, by far the most extensively 
monitored object. 

There is a well-defined radius that characterizes the BLR. A source 
of intrinsic uncertainty involves the likelihood that the BLR is not a 
thin shell of gas at the estimated tblr but a rather thick and stratified 
structure. This means that the power spectrum of spatial frequencies 
in the cross correlation function of continuum and line data can be 
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complex. 

6. The line response is linear. The responsivity of the Balmer lines are gen- 
erally anticorrelated with the incident photon flux (IPronik and Chuvaevl . 



19721 ). Thus, the responsivity varies with distance in the BLR for fixed 



contin uum luminosity and it ca n change with time as the continuum 
varies ( IKorista and Goadl . 120041 ) . This effect might help isolate a par- 
ticular region that is truly virial. 
7. Line emission is roughly isotropic. This assumption is unlikely to hold 
for Balmer lines. In this case the angular patter n of emissivity shoul d 
favor the irradiated face of the line emitting cloud ( iFerland et al.l . Il992l ) . 



More technical problems involve unevenly sampled data, errors in flux cal- 
ibration, normalization of s pectra (different setup, aperture effects), and d i- 



lution by st ellar continuum ( [Peterson! . Il993t iHorne et al.l . I2004J : iBentz et al. 



20061 . 120091 ) . This last problem has been found to ser iously affect e stimates 
of the index in the tblr- L relation. According to iKrolik! ( 1200 ll ). effects 
of a broad radial emissivity distribution (item 3), an unknown degree of 
anisotropy (item 7), and poor sampling in line and continuum monitoring 
can cause additional systematic errors in Mbh determination as large as a 
factor of 3 or more in either direction. 

With reference to item 4, things appear to be different on a timescale 
longer than a couple of years. Looking in more detail at the behavior of 
NGC 5548, and separating the 13 years of data points of continuum and H/3 
flux on a year-by-year basis reveals a distinc t correlation between time delay 
and continuum flux (IPeterson et al.l . |2002| ). The BLR appears to breath 



in response to continuum changes involving a significant change in tblr by 
a factor 3 if we exclude some outlying data points. Even if the line width 
decreases with increasing t B lr (as expected for virial motions), an application 
of Eq. [3] yields a somewhat different mass at each epoch when the continuum 
luminosity is d ifferent - in principle an absurd result. Using the data for H/3 



of NGC 5548 (IPeterson et al.l . 120041 ). we find that Mbh changes by a factor 



4 if the year-by-year data are separately used. If they are averaged together, 
Mbh ~ (1.3±0.3)-10 7 M Q , where the uncertainty is at la confidence level and 
could be taken as indicative of the maximum statistical accuracy achievable 
on the virial product tblrO" 2 using the reverberation method. 

We remark that this degree of accuracy is not possible for the wide ma- 
jority of the sources with tblr determination from reverberation data. A 
major limitation of the reverberation technique involves the need for many 
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spectra of an individual source and spectra spaced across a large enough time 
interval to capture well-defined continuum fluctuations and subsequent line 
responses. This requirement led to a worldwide campaigns tha t resulted in 



reverberation measures for H/3 of more than 60 low z sources ( IBentz et al. 



20091 . 1201CH : iDenney et all |2010| ) . Only a few sources were observed in cam- 
paigns that allowed two or more independent determinations of Mbh- Dif- 
ferences in deriv ed Mbh can be as larg e, a factor 3 for Mark 817 and a factor 



2 for Mark 110 flPeterson et all . 12004 ). 



There is a tight correlation between X-ray variability amplitude and Mbh 
estimates for sources with H/3 reve rberation data. Th e intrinsic dispersion in 
a best-fit linear relation is 0.2 dex ( jZhou et all 120101 ). The method based on 
X-ray variability relies on the assu mption that the X-ray va r iability timescale 



corre lates with source size (e.g., lHayashida et all Il998t iNikolajuk et al 



2004 ). The X-ray method has the advantage that it does not rely on the 
virial assumption. The agreement is very good and the rms scatter quoted 
above is probably representative of a typical statistical accuracy of M B h values 
derived from the reverberation method. 

A few reverberation measures of CivA1549 also exist using IUE (e.g 



1998 



1991 



Gas keil and Sparkdll986l ; iKoratkar and Gaskel]||l991t iTurler and Courvoisier 



and the early Intern ational AG N Watch campaigns e .g., IClavel et al. 



Korista et al.l Il995l ) and HST (IPeterson et all 120041 ) . However, the 



CivA1549 line profile is probabl y composed of one o r more components whose 
broadening is eithe r non- virial (ISulentic et al.ll2007t see also § 17. 2p or affected 
by scattering (e.g., Gaskell and Goosmann . 20081 ). The interpretation of the 
CivA1549 reverberation data is therefore expected to be even more complex 
than the one of H/3. Unambiguous results may be possible only if velocity 
resolved reverberation mapping is carried out. 



5. A more widely applicable technique: -Mbh scaling with luminos- 
ity 

5.1. The tblr - L correlation 

Clearly we need a technique that can provide reasonable estimates for 
the BLR radius in many hundreds or even thousands of quasars. Virial 
estimation appears to be the obvious and only way to accomplish this goal 
but we require a faster way to obtain estimates of tblr for large numbers of 
sources. Reverberation radii are rare and extremely rare for high luminosity 
(logL ^ > 45) quasars which are the sources for which we are currently 
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able to carry out only single-epoch FWHM H/3 measurements beyond z 
0.7. We proceed to correlate all rever beration valu e s with measures of source 
lumin osity. The resultant relation (IKaspi et al.l . |2000| . 120051 ; iBentz et al 



20091 ) provides a secondary estimate of tblr- We assume that the BLR 
scales with luminosity beyond logL « 45 so that the relation defined for 
low luminosity quasars can be directly employed for quasars up to 2 dex 
higher luminosity. It is now standard practice to estimate the BLR radius 
by assuming 



tblr oc (AL; 



(5) 



The exponent in the correlation between so urce luminosi ty L and tblr is 
0.5 £j a £j 0.3, where a ~ 0.7 was obtained by iDibail (Il977l ). The exponent 
has usually been assum ed to be a= 0.5 - 0.7, with a Ri0.52 now considered 



the most reliable value (IBentz et all 120091 ) . Any deviation from this value 
has quantitative effects that are rather modest. If a restriction is made to 
the most likely range, 0.5 £ a 0.6, the effect on M B h estimates is ~ 0.3 
dex over three orders of magnitude in luminosity. 

If we employ H/3, L\ is the specific luminosity at 5100 A in units of erg 
s -iA -1 ; the most homogenous sample available to-date yields the relation 
shown in Fig. HI 

It is important to stress that the relation is mainly based on low z ( £j 0.4) 
quasars and that the high (and low) luminosity ends of the correlation are 
poorly sampled. Con sidering only sour ces in the range 43 £ logL/L Q £ 45 
(i.e. where sources in Bentz et al. 20091 show uniform sampling), we still have 
two orders of magnitude to go in order to reach the most luminous quasars. 
Efforts are now underway for long term monitoring of high- luminosity quasars , 



although they ha v e yet to produce si gnificant results ( Botti et al. . 2010l : 



Kaspi et al.. 


2007; 


Trevese et al.. 


2007) 



Similarly, there is a program at- 
tempting to cover the lower and middle lu minosity range of the correlation, 



the Lick AGN Monitoring Project (LAMP) f lGreene et al.l . l2010al : lBentz et al. 



20101 ; iDenney et al.l . |2010| ). As it is evident from Fig. H] where yearly deriva- 
tions of tblr for NGC 5548 are shown in blue, there will be some irreducible 
intrinsic scatter. Using the L-r BLR relation yields a loss of accuracy for M B h 
estimation; from Fig. H]one infers that the scatter is at least a factor 2 at a 2a 
confidence level at logL m 43.5, although a factor 6 seems possible at lower 
luminosity where data points a may also begin to deviate from linearity. 
Fig. [5] shows FWHM H/3 as a function of luminosity for several samples, 
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including the high- luminosity sample described in iMarziani et al.l ( 120091 ). 
The minimum FWHM H/3 apparently increases with luminosity. The plotted 
lines refer to an tblr- L relation with indices a = 0.52 and a = 0.67, assume 
virial motion and sources at each L radiating at Eddington limit (i.e., the 
limiting curve are not an effect of any flux-limit condition). In principle, 
the trend observed can be used to constrain the exponent of the tblr - L 
co rrelation if larg e samp les of high-L quasars are used. The data presented 
in IMarziani et al.l ( 120091 ) for intermediate z quasars apparently favor a value 
larger than a = 0.52, although the lowest-FWHM points could be due, for 
example, to orientation effects ( §3.3p . 

5.2. Mbh - L relations 

The correlation between tblr and luminosity has been redefined for the 
UV continua luminosity appropriate for using the width of UV lines as a 
virial broadening estimator. This simply means to compute a best fit between 



H/3-based r R LR (e.g., 


McLure and Jarvisl. 2 


002; 


McLure and Dunlopl. 


2004: 


Vestereaard and Peterson!. 


20061: IShen et al. 


. 2008: 


Rafiee and Hall. 


2011). 




Vestergaard and Peterson! (12006) derived mass scaling relations from the 



empirical telr- L relation that have been used to compute Mrh for large 



sampl es of quasars with single-epoch H/3 or CivA1549 observations. IVestergaard and Peterson 
(120061 ) then applied the scaling relations to the objects whose t B lr is avail- 
able from reverberation mapping. They found that the statistical accuracy 
of Mbh mass estimates is low, « ±0.66 and ~ ±0.56 at a ±2cx confidence 
level for Hfl and CivA1549 respectively. Improvements could come from a 
more refined consideration of the virial broadening estimator for single epoch 
observations 



6. Virial broadening and distinguishing quasar populations 

Single epoch measures of FWHM are not necessarily the best estimators of 
the BLR velocity dispersion. The highest precision measure of the virial prod- 
uct involves emission-line width obtained using the cross-correlation function 
centroid (as opposed to the cross-correlation function peak) for the time de- 
lay and line velocity dispersion a (as opposed to FWHM) for the line width. 
Obviously one wants to measure the line width in the variable part of the 
profile. This is the part most likely to be optically-thic k to the Lyman c pntin- 



uum and therefore responding to continuum changes. iPeterson et al.l (120041 ) 



17 



OMarziani et al. 2003 

• Zamfir et al. 2010 

• Marziani et al. 2009 
A Zhou et al. 2006 




J bol 



Figure 5: FWHM as a function of bolometric lum inosity reported for tw o low-z samples 



and the high-luminosity, intermediate z sample of Marziani et al. ( 2009f ). The last sam 



pie is based on h i gh s/ n observations of Hf3 with the IR spectrometer ISAAC at VLT 
( Moorwood et al ., 1998). The shaded area emphasizes an apparent avoidance zone where 
no source is expected if quasars do not radiate super-Eddington, and the exponent of the 
tblr-L scale relation is a w 0.67. Given that the relevant trend is for the minimum FWHM 
as a function of L boi, we show the narrowest sources of a large SDSS sample of NLSyls 
( Zho u et al 1. l2006h . Th ese sources also follow the expected minimum FWHM dependence 
on Lboi ■ Adapted from Marziani et al. ( 20091 ) . 
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Figure 6: H/3 profiles in spectral types A2 and Bl that are typical of Pop. A and B 
respectively. The displayed median spectra have been computed from a large sample of 
SDSS high s/n, low-z quasars (Sulentic et al., in preparation). 



found that that the rms a minimizes the random component of errors in 
reverberation-based Mbh measurements (with an absolute minimum around 
30% for NGC 5548, as discussed in §1]), so that rms a is probably the best 
choice. 

Which measurement is better in general for single epoch spectra, FWHM 
or cr? FWHM is very straightforward; however it is sensitive to substructures 
in the profile (e.g., double peaked profiles), and presence of inflections that 
can yield unstable values depending on continuum placement. On the con- 
verse, a has the significant disadvantage of diverging in profiles with promi- 
nent line wings (and to be — > oo for Lorentzian profiles). Especially in the 
case of low s/n, a values are affected by large errors. 

On e must at least co nsider that the ratio FWHM/a changes with line 
width (ICollin et all 120061) and avoid mixing all sources together in large sam- 



ples (e.g., iNagao et al.l . 120061 ) . The change occurs when radia tive and gravi- 



tation al forces appear to balance (Eq. |4j: L/L Edd ^ 0.15±0.05 ( iMarziani et al. 



2003bf) for iV c = 10 23 cm 2 , probably an appropriate average value of column 
density. 
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6.1. H/3 in Population A Sources 

Pop. A profiles (FWHM H/3< 4000 km s" 
nent Lorentz function (Gaussians yield larger resid uals; IVeron-Cetty et al 



are well fi t by a single compo 



200ll ; ISulentic et al.l . |2002| ; iMarziani et all l2003bl ) and tend to show lower 



amplitude variations and less change in profile shape when continuum varie s 
(compared to Pop. B(roader) sources; iKlimek et al.l 12004 lAi et al.ll2010h . 
One can describe them as Narrow Line Seyfert 1-like sources. The biggest 
challenge for extracting an accurate FWHM measures from Pop. A profiles 
involves careful subtraction of the usually strong Fell emission that afflicts 
and broadens the red side of the profile. Line profile shifts or asymmetries 
are usually small. 

6.2. H/3 in Population B Sources 

The typical Pop. B H/3 line (FWHM H/3 > 4000 km s" 1 ) is more com- 
plex and requires two Gaussian component models to adequately describe 
the profile (Fig. [6]). The narrower, broad component is likely the most re- 
liable virial estimator while the VBC component showing a redshift of 1000 
- 2000 km s" 1 and FWHM 10000 km s~ x is unlikely to be virial. The nar- 
rower component then is viewed as the classical BLR component which is 
the only one seen in Pop. A sources. As for Pop. A, usually shifts are small 
although one must consider that a fraction o f Pop. B s ources show shifts 



that are significant relative to the line width; IZamfir et al.l 120101) . It appar- 
ently reverberates more strongly than the VBC (e.g., ISulentic et al.l . |200( 
Korista and Goadl . 12004 ) . The VBC can affect both estimates of tblr and 



virial velocity for almost half of quasars. It must be taken into account dur- 
ing reverberation studies or it will slow down and/or blur the CCF results; 
it must be considered in measures of FWHM H/3 or the virial velocity can 
be seriously overestimated (factors of 2 or 3). The relative strengths of the 



two H/3 components appear s to depend on source luminosity (IMarziani et al. 



20091 ; iMeadows et all 1201 ll ) with the VBC becoming more dominant in high 



luminosity quasars. This is obviously important for studies of r blr and M B h 
as a function of L and z. As spectral s/n decreases the two components are 
unlikely to be recognizable with the result that FWHM H/3 will be system- 
atically overestimated. 
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7. Mass estimates for large samples 



We must examine all strong broad lines as potential virial estimators. 
In the following, this review focuses on attempts to estimate Mbh for large 
numbers of quasars over as wide as possible range of redshift and source 
luminosity. Opinions will vary about how far we have succeeded in this 
quest. 

Whether Pop. A or B, the number of sources with reliable FWHM H/3 
measures becomes increasingly small above z ~ 0.7. One is left with three 
options if one wishes to study sources at higher redshift: (1) follow H/3 into 
the infrared, or (2) adopt other broad lines (e.g. MgllA2800 or CivA1549) as 
surrogate virial estimators. If one wishes to study lower luminosity sources 
at low redshift then option (3) involves using databases like SDSS where 
cont extual binning of many source spec t ra ca n yield high s/n profiles (see 
e.g. IZamfir et al.l . |2010| ; iMeadows et al.l . 1201 ll ). The availability of infrared 
spectrographs at 4m+ class telescopes makes option (1) possible for tens or 
even hundreds of sources. By "possible" we mean that it becomes possible to 
obtain spectra with s/n (~ 20 in continuum near H/3) comparable to lower 
redshift optical measures for the brightest individual SDSS (or other) spectra. 
Using ESO VLT and IR spectrometers like ISAAC one can presently obtain 
such spectra for sources at the bright tip of the quasar luminosity function 
up to z ^ 3.8 with exposure times of ~40 minutes. The availability of IR 
windows determines the source redshift distribution using H/3 but reasonable 
coverage over the range z — 1 — 3.8 is possible. The advantage of this 
approach is uniformity, same reduction procedure and same virial estimator. 

The SDSS database is a seductive source of data for more than 100000 
quasars offering the possibility to derive instant Mbh estimates. Beyond the 
brightest « 500 sources this is an illusion as the decline in s/n makes reliable 
FWHM estimates increasingly difficult. If all quasar spectra were basically 
the same then one could proceed by fitting profile models, determined for 
the brightest sources, to large numbers of noisy spectra. We now know that 
spectra of quasars are less similar than are stellar spectra along the main 
sequence of the H-R Diagram. Dissimilarities include: 

1. internal broad (and narrow) line shifts, 

2. the line profiles at a fixed redshift or source luminosity show large 
differences as detailed earlier, 

3. Fell contamination varies widely and directly affects the FWHM mea- 
sure. 
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Quite aside from arguments about the validity of any specific line as a virial 
estimator, the SDSS database of spectra (all but 100s) can be used but only if 
the data are bi nned into some well d efined context (for example the spectral 
type defined bvlSulentic et al.ll2002l ). or some innovative approach is applied 



e.g., iRafiee and Halll . 1201 if ). 



1.1. Mass estimation using Mgll\2800 



McLure and Dunlopl ( 120041 ) computed virial black hole mass estimates for 



a sample of more than 10000 quasars in th e redshift interval 0.1 < z < 2.1 
drawn from the SDSS Quasar Catalogue of ISchneider et al.l (120031 ). The H/3 
line was used up to z ~ 0.7, and MgllA2800 beyond. The change in virial 
estimator did not introduce any obvious discontinuity. Also no evidence 
of large masses exceeding the Mbh of locally detected dormant black holes 
were found. This behavior is not surprising since H/5 and MgllA2800 are both 
LILs thought to be produced mostly in the same emitting region. There is 
evidence that t he MgllA2800 FWH M distribution might be narrower than 



the one of H/3 (ILabita et al.l . |2009|), and that individual MgllA2 800 FWHM 



values are systematically narrower than H/3 (jWang et al.l . 120091 ). If so (and 



if the presence of narrow absorptions on the blue side of MgllA2800 are not 
creating a spurious effect), MgllA2800 might be more appropriate because its 
single epoch profile may be a better approximation of the rev erberating part 



of H/ 3 isolated on the rms spectrum, as also suggested for Fell (jSulentic et al. 
20061 ). 



1.2. (Dangerous) mass estimation using Cw\1549 

Much more controversial has been the use of CivA1549 as a virial broad- 



ening estimator (IBaskin and Laorl . 120051 ; iNetzer et al.l . 120071 ) 



The most serious effect is related to the str ong blueshift with respect to 



the re st frame frequently observed in CivA1549 ( lGaskeHl . ll982tlRichards et al. 



201ll ). This sugges ts that either t he high-ionizatiqn lin e emitting gas is out- 
flowing in a wind ( jGaskelll . Il982l ; iMarziani et al.l . 119961 ) in at least a signifi- 
cant fraction of quasars and hence not in virial equilibrium, or that there is 
some other broadening mechanism , such as electron scattering or Rayleigh 
scattering ( IGaskell and Goosmannl . 120081 ). In either case CivA1549 is not a 
reliable mass indicator. This fraction may even increase with redshift and 
lumino sity, judging from th e observed CivA1549 profile in a high z quasar 
sample iBarthel et al.l ( 119901 ) . The effect is very pronounced if we consider a 
prototypical source believed to accrete close to the Eddingon limit, I Zw 1 
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1500 1550 1600 1 1.5 2 2.5 

Rest frame wavelength [A] Log(W(CIVX1549 BC )) [A] 



Figure 7: Left: CivA1549 and H/3 profile of the prototype Pop. A source I Zw 1. The 
grey line represents the original profile, the black one the scaled H/3 profile and the blue 
solid line the residual after subtraction of the scaled H/3 profile. The ratio between A/bh 
computed using FWHM CivA1549 and FWHM H/3 as virial broadening estimators as 
a function of W(CivA1549). Filled circles: Pop. A radio-quiet sources; filled squares: 
Pop. A radio- quiet; open circles: P op. B radio-quiet; open squares: Pop. B radio-loud. 
Adapted from lSulentic et all (|2007l) . 
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(Fig. [7]). The CivA1549 profile can be modeled as an the unshifted almost 
symmetric H/3 profile plus a prominent blueshifted component. The broad- 
ening, in this case, is dominated by a shifted non-virial component. Since 
M-qk oc (Sv ) 2 , the effect can reach a factor 10 on Mbh, as shown in the right 
panel of Fig. For Pop. B sources, there is no obvious systematic ef- 
fect, although the Mbh derived from CivA1549 is sometime much lower than 
the one derived from H/3, probably because of significant CivA1549 narrow 
emission (CivA1549 N c)- 

The very existence of CivA1549nc has been a contentious issue. Radio- 
loud sources often show prominent CivA1549 cores o f width ~ 2000 km s _1 , 
but such emission has been as cribed to CivA1549rc (Wills et all 119931 ) and 
to an intermediate line region (IBrotherton et all . Il994l ) . We have interpreted 
this feature as part of the NLR, and recent detection of na rrow CivA1549 
in high z type-2 sourc e s sup ports this interpretation (e. g., iMainieri et al. 



20051 ; ISevergnini et all 120061 ) . Provided that a smooth density gradient exist 



(i.e., that there is a significant amount of gas at n e ~10 6 cm -3 , we ex- 
pect CivA1549 emission to be strong (CivA1549 emissivity is oc n 2 while 
[Olll] AA4959,5007 should become collisionally quenched). If velocity disper- 
sion increases with decreasing distance from the central continuum source, 
as expected if there is at least rough virial equilibrium, then it is some- 
what natural to expect FWHM(CiyA1549 Nn ) > FWHM£[Qm] AA4959,5007) 



(IDe Robertis and Osterbrockl . Il986l ; lAppenzeller and Oestreicherl . Il988l ; ISulentic and Marziani 



19991). 



A comparison of F WHM H/3rc and CiyA1549 line shift measur es with 



those of other workers (IWarner et all . l2004j ; iBaskin and Laorl . 120051 ) shows 
that significant and systematic differences exist between the measures. FWHM 
CivA1549bc measures begin to deviate towards smaller values at FWHM 
CivA1549 ~ 4 000 km s -1 where w e believe that the NLR becomes impor- 
tant (Pop. B; iBachev et~ai1 l2004h . However, considering the strong Bald- 
win effect that is observed in [Olll]A5007 (mainly du e to the increasing 
rarity of [Oin]A5007 -strong objects at high luminosity iNetzer et all 12004 ; 



Kovacevic et alll2010l ; Sulentic et al., in preparation), the relevance of CivA1549 



NC 



could be small in high z quasars, for which CivA1549 becomes the only line 
available for estimating the virial broadening from optical observations. 

Concluding, the FWHM of CivA1549 is a poor virial broadening estima- 
tor subject to systematic biases. Alternatively the intermediate ionization 
lines (especially Sim] A1892 and AI111AI86O) could in principle be used even 
if blended. Intermediate-ionization lines like Sim] A1892 and AI111AI86O show 
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negligible blueshifted, non virial emission, we ak VBC and a profile roughly 
consistent with the broad component of H/3 (IMarziani et al.l . 120101 and ref- 
erences therein). 



8. Photoionization Methods 



Historicall y, the f i rst at tempt to derive Mbh from photoionization condi- 
tions is due to iDibail ( 119771 ). as mentioned in ^TJ The underlying assumption 
was of course that all AGN were in similar physical conditions, so that the lu- 
minosity of H/3 could be related to the BLR size assuming an identical value 
of the line emissivity. The application of Eq. [I] needs an estimate of the fill- 
ing factor. If spherical symmetry is assumed, /f can be ea sily computed from 
the c overing; factor f c that is observationally constrained (IWandel and Yahill . 
1985h . 

The physical conditions of photoionized gas can be described by hydrogen 
numeric density or electron density, hydrogen column density iV c metal- 
licity, shape of the ionizing continuum, and the ionization parameter U. The 
latter represents the dimensionless ratio of the number of ionizing photons 
and the total hydrogen density, ionized and neutral. If we know the product 
of rift (or electron density n e ) and U, we can estimate tblr from Eq. |2j U 
can be written in terms of tblr yielding an estimate of the BLR independent 
of direct reverberation estimates: 



r BLR 



f ¥dv 
AmJJriftC 



oc 



(Un H ) 



tQ(hV 



(6) 



The dependence of U on r blr was used by lPadovani and Rafanellil (119881 ) 
to compute centra l black ho l e mas ses assuming a plausible average value of 
the product nJJ . iPadovanil ( 119891 ) derived an average value U ■ n e « 10 9 8 
from several sources for which tblr was available from reverberation map- 
ping and for which the number of ionizing photons could be measured from 
multiwavelength observations. The average n e U value was then used to com- 
pute black hol e masses for a much larger sample of Seyfert 1 galax ies and 
low- z quasars (IPadovani and Rafanellil . Il988l ; iPadovani et al.l . Il990l ). Mul- 
tifrequency data were used to define the shape of the ionizing continuum 
for each individual source. IWandel et al.l (119991 ) compared photoionization 
method results with those obtained via reverberation mapping and found a 
very good correlation between the two mass estimates. 
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Recently, the photoio nization method of iDibai (11984T) based on H/3 lumi- 
nosity was reconsidered (IBochkarev and Gaskelll . 120091 ). The original Dibai 
data yield very good agreement with Mbh estimates from reverberation- 
mapped sources. 

A different photoionization method could be applicable to high redshift 
quasars. Emission lines originating from forbidden or semi-forbidden tran- 
sitions become collisionally quenched above the critical density and, hence, 
weaker than lines for which collisional effects are still negligible. The Aim A 1860 
/Sim] A1892 ratio is well suited for exploring the density range 10 11 — 10 13 
cm -3 which corresponds to the densest, low ionization emitting regions likely 
associated with the production of Fell. The ratios Sill] A1814 /Sim] A1892 and 
SilvA1397/ Sim] A1892 are independent of metallicity and sensitive to ion- 
ization. Conversely the ratio ClV A1549/Siiv A1397 is mainly sensitive to 
metallicity. 

We computed a m ultidimensional grid of CLOUDY ( IFerland et al.Ul998[ ) 



simulations (see also iKorista et al.l . 119971 ) in order to derive U and nn- Com- 
putation of constant value contours in the theoretical U vs. riu plane of 
CLOUDY simulations shows convergence towards a low ionization plus high 
density range. Application of this method to the NLSyl I Zw 1 and NLSyl- 
like sources at high redshift provi des not only an indep endent estimate of nil 
but also of n H , U, and metallicity ( jNegrete et al.l . 120101 ) . r blr values obtained 
with the photoionization method can be compared to the ones obtained with 
sources reverberation mapping. The two sets of values agree better than the 
values estimated using the tblr- L correlation - and th i s emp loying the very 
sample used to define the tblr- L correlation (jNegretd . 1201 ll ). 

The large uncertainty associated with the scaling relation that employs 
CivA1549 (present even if the CivA1549 virial component can be properly 
isolated) makes preferable a one-by-one Mbh determination based on physical 
properties of an emitting region that remains self-similar. The particular 
value of the photoionization method is that it can be applied over the full 
redshift range where quasars are observed. 
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9. Black Hole Mass and Host Galaxy Mass 



The correlation of nuclear black hole mass with stellar bulge velocity 



dispe rsion in nearby galaxies ( iFerrarese and Merrittl . I200CH ; iGebhardt et al 



20001 ) suggests a constant M B H/M bulge ratio. This result has been extensively 



used (or should we better say abused?) to derive quasar M B h from the ratio 
M Rn/Mhnigfi or from a pr o xy of the stellar velocity dispersion. 



Ferrarese and Merritt (2000) relied on maser optical emiss i on lin e mea- 



surements. IGebhardt et al.l ( 120001 ) and iMerritt and Ferraresd (120011) added 
mass e s obtained from the so- called stellar dynamics method (e.g. iKormendyl . 



19931 ; iMagorrian et al.l . Il998l ). Broadly speaking, these Mbh determinations 



have the advantage that the line emitting region is partly resolved, so that 
a detailed analysis of geometry and orientation effects can be carried out. 
It is still debated whether stellar dyn amics Mbh determin ations might in- 



troduce a bias in the measured Mbh (IGultekin et al.l . |2009| ). However, it is 



more relevant in the the quasar context to assess whether a significant Mbh~ 
Mb u ig e correlation exists at all, or whether it results from a bias due to the 
difficulty of detecting a black h ole whose spher e of in fluence is much smaller 
than the telescope resolution (IGultekin et al.l . 1201 ll ). In other words, the 
correlation may refer only to the maximum Mbh for each Mb u i ge . The most 
likely condition is that the bias is mass-dependent, and most severe at low 
Mbh (i-e., higher mass black holes are more easily detected). The true (unbi- 
ased) aspect of the relation could be tight at large Mbh but with increasing 
dispersion toward lower masses. 

An additional complication is that the relation bet ween Mbh and a* ap 



pears to be different for barred and non-barred hosts ( iGraham et al.l . 12011 



and references therein). Intrinsic dispersion at low Mbh and dependence on 
morphology could overshadow several conflicting claims e.g ., a nonline ar rela- 
tion between the Mbh and the bulge mass for PG quasars ( iLaorl l200ll ): Mbh 
oc M^ u ^ °' 14 . Another ca s e is th e claim that NLSyls, often host in dwarfish 
galaxies ferongold et al.l . 1200 ll ) and barred spirals (ICrenshaw et al.l. 12 003 



Ohta et al.l. 120071 ) poss ess under-massive black holes (e.g., iMathur et al 



20011 ; IChao et all . |2008|). 



A similar difficulty migh t apply also to the correlat i on be tween black hole 
mass and bulge luminosity ( IKormendy and Richstond . Il995l): the dispersion 
of da t a at low-Lhiiw c ould make any correlation ill-defined (IGaskell and Kormendyl . 
20091 ; iGaskelll . l2010al ). The L hu \ ge - M B h relation suffers from the additional 



problem that Seyfert galaxy hosts are brighter than normal galaxies for a 
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given value o f their velocity disp ersion, perhaps as a result of younger stellar 
populations (INelson et al.l.l2004j), even if the use of IR luminosity can greatly 
reduce the scatter (IMarconi and Huntl . 120031 ) . 

Radio-quiet quasars seem to follow the established Mbh-c* re lation up to 
z ~0 . 5, with a modest evoluti on in the redshift range 0.5 ^ z ^ 1 ( [Shields et al. 



2003 



Salviander et all 120071 ). However, the M B u/ M hu \ ge ratio estimated at 
low z should not be mistaken for an universal constant. Studies at high red- 
sh ift find Mrh and Mh „]p- e values that indicate overmassive black holes at high 
z ( iTargett et al.l . l201ll . and references therein). Evolution in the Mbk/ M^ge 
ratio casts doubts on the physical implications of the relation at \ow-z. 



10. Using [Olll]A5007 as a proxy for stellar velocity dispersion 

If the M-Qn/M^ige ratio can be derived with good accuracy, then M B h 
could be computed using a proxy for a*. There are two major difficulties 
with using FWHM [Olll]A5007 as a proxy f or o*. While a M BH - FWHM 
[Oin]A5007 correlation exists (INelsonl . |2000| ). the scatter is large. The de- 
rived v alues can be considerab ly higher than those calculated using FWHM 
H/3bc ( IMarziani et al.l . l2003b[ ). On the one hand, the narrowest profiles of 



[Olll]A5007 appear to be sub-virial in the gravitational potential of the host 
bulge, probably because the line emitting gas is constrained in non- spherical 
geometries (see the discussion in iGaskelll l2009al and references therein) . On 
the other hand, low-W([Oin]A5007) sources can show FWHM [Olll]A5007 
FWHM(H/3), seriously challenging the assumption of gravitational dom- 
inance by the bulge potential. Bl ueshifted [Qui] A5007 profiles are observed 
and likely arise in outflowing gas (jZamanov et al.l . 120021 ). possibly associated 
with a disc wind. The distribution of [Olll]A5007 velocity shifts with re- 
spect to [Oll ]A3727 indicate s that the majority of sources are significantly 
blue-shifted ( IHu et al.l . 12008 ). The narrow line region (NLR) in blueshifted 
sources is not likely to be dynamically related to the host-galaxy stellar bulge. 
This is even more true for a minority of sources w ith At> r ^ — 300 km s _1 
(the so-called blue-outliers of IZamanov et al.ll2002l ) whose NLR may be very 
compact. This points to a limiting W([Olll] A5007) ~ 20 Abelow which 
FWHMQOlll] A5007) emission probably ceases to be a useful mass estima- 
tor. Only large W([Olll] A5007) radio-quiet Pop. B sources may have very 
extended NLRs whose motions are dominated by the stellar bulge. Very ap- 
propriately, the width of the [Sll]AA6716,5631 or [Oll]A3727 doublet has been 
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proposed as an alternative t o [Olll]AA4959,5007 width (IKomossa and Xul . 
20071 : ISalviander etaD . l2007h . 



11. Masses from Double Peakers 



A small fraction of AGN exhibit ver y broad and double peaked LILs. The 

Ha BC profile is strik i ngly p eculiar (see Sulentic et al.llioOOa ; Eracleous and Halpern 
20031 ; IStrateva et al.l 120031 for examples). Prototype double-peakers Arp 
102B, 3C 390.3, and 3C 332 have been monit ored for more than 20 years 
(ILewis and Eracleousl . 120061 ; iGezari et al.l . 120071 ) . A common property of dou- 
ble peakers involves variability of the profile shape on timescales of months 
to years. This slow systematic variability of the line profile shows similar 



timescales to dynamical cha nges expected in an accretion disk (lEracleous and Hafpernl . 



2003 



Newman et al.l . Il997l ). Five year monitoring of H/3 (IShapovalova et al 



20011 ) supports rejection of the binary BH hypothesis (Gaskell, 1983 ) fo r 
3C390.3, on the basis of the masses required ( IHalpern and Filippenkol . Il988l ). 
Hot spots, spiral waves and elliptical accretion disks have als o been suggested 
( IStorchi-Bergmann et al.l . 120031 ; ILewis and Eracleousl . 120061 ) . 

If one can detect changes in the line profile induced by the hot spot, it 
is in principle possible to derive the period of the hot spot and hence to 
estimate tblr and Mbh in physical units. In the case of Arp 102B one ob- 
i (2.2 ± 0.7) • 10 8 M for the black hole, and a distance of 
500 R g ) for the hot spot ( iNewman et all 119971 ). Derivation 



tains a mass 
4.8-10" 3 pc (f 

of tblr in physical units removes the degeneracy introduced by the Keple- 
rian velocity field (i.e., the velocity scales as (Mbh / r)~ 05 , and disk model 
profiles yield a distance normalized to R g ). Unfortunately, gravitational and 
transverse redshifts show the same dependence oc(Mbh / r)~ 5 at first or- 
der. A gravitational redshift will not help us solve for Mbh and r in physical 
units; detection of a radial velocity change with time will do it. Monitoring 
has been pursued for sources showing single-peaked line profiles under the 
assumption t hat the disk contribution is masked by emission from a wind 



or hot spots ( ]Bon et al.l . 120091 ; IJovanovic et al.l . 120101 ) . This approach pro 



vided a new estimate of the black hole mass for NGC 4151 (Bon et al., in 
preparation). 



12. A binary black hole revival 

A binary black hole model was originally prop osed to e xplain sources 
showing double-peaked line profiles like Arp 102B ( iGaskelll . Il983l ). Lower 
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Rest frame 



Rest frame 




Figure 8: Sketch illustrating 2 examples of most- frequently observed H/3 profiles in Pop. B 
sources (upper half) and the profiles that might be signatures of bina ry black holes. The 
SDSS 153636.221+044127.0 U/3 line profile (|Boroson and Lauerl liooih corresponds to the 
lower-left case. 



limits for plausible orbital periods based on the absence of peak radial velocity 
changes would req uire supermassive binary black holes with total masses i n 



excess of 10 10 M ( iHalpern and Filippenkol . Il988l : IShapovalova et al.L 1200 ll ) 



Such large Mbh values are difficult t o reconcile with the m aximum Mbh 



19971 ) and with the 



values observed in the local Universe (lEracleous et al 
Mbh - bulge mass relation. 

The search for supermassive binaries with sub-parsec separat ion has un- 
dergone a revival (lEracleous et al.L 1201 ll ; iTsalmantza et al.L 120111 ) since such 
binary black holes are expected in scenarios involving bla ck hole and bulge 



coevolution ( jBegelman et all Il980t IVolonteri et al.L 120091 ) . A candidate in 



volves the quasar SDSS J153636. 221+044127.0 that shows line features pos- 
sibly indicating two broad-line systems separated in radial velocity by 3500 
km s _1 , and a system of unresolve d absorption lines at intermediate radial 
velocity ( Boroson and Lauerl. 120091) . A disk model with off- axis illumination 



has also been proposed (IGaskelll . l2010bl id) for this source. 
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It is difficult to use a single source with almost unique^ line properties 
to argue for a phenomenon (i.e. binary BH) expe cted to be common among 
quasars. Peculiar shifts of the kind identified by iBoroson and Lauerl (120091 ) 
are restricted to Pop. B sources and the most frequently observed profile 
showing a redward asymmetry and a peak displacement to the red is clearly 
unsuitable as a binary BLR signature where the NLR indicates the rest frame 
(Fig. [S]). If we interpret the H/3 BC and VBC as two different BLRs asso- 
ciated with the binary, and we assume that the line emitting gas forms two 
bound systems, then the largest shifts should be observed in the BC, the con- 
trary of what is most frequently observed (the VBC is redshifted by ^ 1000 
km s _1 ). Other systems have been selected as possible candidates, notably 
the ones showing oppo site peak and line base displacement in radial velocity 
( IShen and Loebl . |2010| ) , but altogether they account for no more than a few 
percent of quasars. 



13. Toward higher and higher redshifts: the evidence for a turnover 

FWHM H/3 measures can provide good Mbh estimates out to z « 3.8. 
8m-class telescopes yield high s/n measures for the high luminosity tail of the 
quasar optical luminosity function. Sampling over a wider part of the optical 
luminosity function (down to Seyfert 1 luminosities (my «-23) becomes pos- 
sible with the next generation of large telescopes. BH mass estimates beyond 
z ~ 3.8 or simply larger samples beyond z ~ 0.7 requires use of H/5 surrogate 
lines. The two most used surrogates are MgllA2800 and CivA1549. CivA1549 
cannot be trusted but MgllA2800 may be able to serve as a surrogate virial 
estimator for the highest redshift quasars currently known. Spectra can be 
obtained for z ~ 6 quasars in the K band ( IKurk et all 120071 ). 

The best Mbh estimates out to z rj 4 show no evidence of a turnover which 
would reflect the epoch when the largest black holes were still growing. In- 
stead we see constant Mbh upper limit near log Mbh ~ 9.7 if we trust in part 



(which we do not) the measurements based on CivA1549 (IShen et a l.. 2008 



Labita et al.l . 120091 ) . There may be a change at higher redshift (ITrakhtenbrot et al. 



201 ll ) if we consider only Mbh values obtai ned from H/3 and M gllA2800 mea- 



sures. Fig. [9] combines a low- z sample ( jZamfir et al.l . 120101 ) with samples 
using H/3 observed with IR spectrometers (to z ~ 3.5) and MgllA2800 us- 
ing optical and IR measures. We see a possible turnover in estimated M B h 



*A twin of SDSS J153636. 221+044127.0 is shown in Fig. 1 of lGaskelll (|l983h . 
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Figure 9: Mbh versus z for 



a low-z sample (grey dots 
samples. Red circles: 



Zamfir et al. 201 



and 



several i mtermediate t o high z samples. Re d circles: | Ma rziani c t al.l <|2009f ) : open 
squares: Dietrich et al. ( 20091 ): open triangles: Shemmer et all (|2004l ): filled pentagons: 
Netzer et al.l (12007 ): filled squares: Trakhtenbrot et al. ( 201l[ ): open starred octagons: 



Willott et all (|2010l ): filled octagons: iKurk et al.l (|2007l ): lar ge spot at z sa 7 



quasar whose discovery was announced in late June 2011 
dashed line marks A/bh = 5-10 9 M . 



the high-z 
(|Mortlock et all . l20lTh - The 
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at the highest redshifts although some care is needed in interpreting these 
results. Most sources in the range 1 £ z £ 2.5 were selected from the bright- 
est quasars in the Hamburg ESO survey and are the most luminous quasars 
known in that redshift range. Observations at very high redshift refer to 
much fainter quasars. We must restrict our attention to the high-end of the 
mass distribution when we evaluate the significance of the turnover. Count- 
ing sources with masses in the ranges 9.25 < log M B h< 9.75 and 8.75 < log 
Mbh< 9.25, we find that the ratio of the numbers of less-massive to more- 
massive sources at redshift ^ 4 is lower than for the samples at 1 z 3.8. 
A simple application of Poisson statistics to these ratios confirms a real trend. 
Given the different sample selection criteria at different redshifts we believe 
that more data are needed before the turnover can be regarded as established. 

As a final consideration we note that the computed Mbh may not be 
critical for concordance cosmology, since black holes can grow to the observed 
masses in a duty cycle that is significantly shorter than the age of the Universe 



at z ~ 6 according to iTrakhtenbrot et al.l (1201 if ) . 



14. Conclusion 

Mbh computation techniques for large samples of quasars are rough and 
the lack of accuracy in Mbh estimates is serious. There are several areas that 
could lead to significant improvement: 

• a significant reduction in scatter could be achieved by more careful 
selection of virial broadening estimators (best are H/3 and Mgll); 

• a second factor is related to know ledge of the BLR structure that is 



still hotly debated ( jGaskelll . l2009bl ). There is evidence that Pop. A and 
B sources show different BLR structure and kinematics. Significantly 
different / values are likely associated with the two populations; 

• photoionization methods should be favored over methods based on the 
tblr-L correlation. 

Considering the large scatter introduced by uncertainties in the factors en- 
tering the virial relation it is still not surprising that Mbh estimates and 
those derived by randomly reassigning the quasar broad-line widths to dif - 



ferent objects show such similarities in the Mbh vs. z plane ( jCrooml . 1201 it ) . 
However this provocative result may not stand for long. 
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The work was presented as an invited talk at special workshop " Spectral 
lines and super-massive black holes" held on June 10, 2011 as a part of 
activity in the frame of COST action 0905 "Black holes in an violent universe" 
and as a part of the 8 th Serbian Conference on Spectral Line Shapes in 
Astrophysics. We are indebted to Martin Gaskell for discussions and many 
insightful suggestions. We also acknowledge with gratitude the hospitality 
and good organization of the Conference in Divcibare (Luka, Dragana, Darko 
and all the others of the organizing committee). 
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